Gene regulatory networks: a coarse-grained, equation-free approach to multiscale computation

We present computer-assisted methods for analyzing stochastic models of gene regulatory networks. The main idea that underlies this equation-free analysis is the design and execution of appropriately-initialized short bursts of stochastic simulations; the results of these are processed to estimate coarse-grained quantities of interest, such as mesoscopic transport coefficients. In particular, using a simple model of a genetic toggle switch, we illustrate the computation of an effective free energy and of a state-dependent effective diffusion coefficient that characterize an unavailable effective Fokker-Planck equation. Additionally we illustrate the linking of equation-free techniques with continuation methods for performing a form of stochastic "bifurcation analysis"; estimation of mean switching times in the case of a bistable switch is also implemented in this equation-free context. The accuracy of our methods is tested by direct comparison with long-time stochastic simulations. This type of equation-free analysis appears to be a promising approach to computing features of the long-time, coarse-grained behavior of certain classes of complex stochastic models of gene regulatory networks, circumventing the need for long Monte Carlo simulations.Comment: 33 pages, submitted to The Journal of Chemical Physic

Biochemical Network Stochastic Simulator (BioNetS): software for stochastic modeling of biochemical networks

BACKGROUND: Intrinsic fluctuations due to the stochastic nature of biochemical reactions can have large effects on the response of biochemical networks. This is particularly true for pathways that involve transcriptional regulation, where generally there are two copies of each gene and the number of messenger RNA (mRNA) molecules can be small. Therefore, there is a need for computational tools for developing and investigating stochastic models of biochemical networks. RESULTS: We have developed the software package Biochemical Network Stochastic Simulator (BioNetS) for efficiently and accurately simulating stochastic models of biochemical networks. BioNetS has a graphical user interface that allows models to be entered in a straightforward manner, and allows the user to specify the type of random variable (discrete or continuous) for each chemical species in the network. The discrete variables are simulated using an efficient implementation of the Gillespie algorithm. For the continuous random variables, BioNetS constructs and numerically solves the appropriate chemical Langevin equations. The software package has been developed to scale efficiently with network size, thereby allowing large systems to be studied. BioNetS runs as a BioSpice agent and can be downloaded from . BioNetS also can be run as a stand alone package. All the required files are accessible from . CONCLUSIONS: We have developed BioNetS to be a reliable tool for studying the stochastic dynamics of large biochemical networks. Important features of BioNetS are its ability to handle hybrid models that consist of both continuous and discrete random variables and its ability to model cell growth and division. We have verified the accuracy and efficiency of the numerical methods by considering several test systems

Mathematical model reveals role of nucleotide signaling in airway surface liquid homeostasis and its dysregulation in cystic fibrosis

The intrapulmonary airways conduct air to the alveoli and are defended from inhaled pathogens by a highly regulated protective system of mucus, cilia, and liquid. In healthy lungs, a well-hydrated mucus layer is cleared by cilia from airway surfaces. In cystic fibrosis (CF), airway surfaces are dehydrated, leading to a failure of cilia-mediated mucus clearance and accumulation of pathogen-infected mucus. In this study, we created a mathematical model of airway surface liquid regulation in normal and CF cells and used the model to investigate a potential therapy to rehydrate CF airways and restore proper mucus clearance

Variable-free exploration of stochastic models: a gene regulatory network example

Finding coarse-grained, low-dimensional descriptions is an important task in the analysis of complex, stochastic models of gene regulatory networks. This task involves (a) identifying observables that best describe the state of these complex systems and (b) characterizing the dynamics of the observables. In a previous paper [13], we assumed that good observables were known a priori, and presented an equation-free approach to approximate coarse-grained quantities (i.e, effective drift and diffusion coefficients) that characterize the long-time behavior of the observables. Here we use diffusion maps [9] to extract appropriate observables ("reduction coordinates") in an automated fashion; these involve the leading eigenvectors of a weighted Laplacian on a graph constructed from network simulation data. We present lifting and restriction procedures for translating between physical variables and these data-based observables. These procedures allow us to perform equation-free coarse-grained, computations characterizing the long-term dynamics through the design and processing of short bursts of stochastic simulation initialized at appropriate values of the data-based observables.Comment: 26 pages, 9 figure

From physics to pharmacology?

Over the last fifty years there has been an explosion of biological data, leading to the realization that to fully explain biological mechanisms it is necessary to interpret them as complex dynamical systems. The first stage of this interpretation is to determine which components (proteins, genes or metabolites) of the system interact. This is usually represented by a graph, or network. The behavior of this network can then be investigated using mathematical modeling. In vivo these biological networks show several remarkable (and seemingly paradoxical) properties including robustness, plasticity and sensitivity. Erroneous behavior of these networks is often associated with disease. Hence understanding the system-level properties can have important implications for the treatment of disease. Systems biology is an organized approach to quantitatively describe and elucidate the behavior of these complex networks. This review focuses on the progress and future challenges of a systems approach to biology

In Silico Generation of Alternative Hypotheses Using Causal Mapping (CMAP)

Previously, we introduced causal mapping (CMAP) as an easy to use systems biology tool for studying the behavior of biological processes that occur at the cellular and molecular level. CMAP is a coarse-grained graphical modeling approach in which the system of interest is modeled as an interaction map between functional elements of the system, in a manner similar to portrayals of signaling pathways commonly used by molecular cell biologists. CMAP describes details of the interactions while maintaining the simplicity of other qualitative methods (e.g., Boolean networks)

Compression and dilation of the membrane-cortex layer generates rapid changes in cell shape

A cyclic process of membrane-cortex compression and dilation generates a traveling wave of cortical actin density that in turn generates oscillations in cell morphology.Rapid changes in cellular morphology require a cell body that is highly flexible yet retains sufficient strength to maintain structural integrity. We present a mechanism that meets both of these requirements. We demonstrate that compression (folding) and subsequent dilation (unfolding) of the coupled plasma membrane–cortex layer generates rapid shape transformations in rounded cells. Two- and three-dimensional live-cell images showed that the cyclic process of membrane-cortex compression and dilation resulted in a traveling wave of cortical actin density. We also demonstrate that the membrane-cortex traveling wave led to amoeboid-like cell migration. The compression–dilation hypothesis offers a mechanism for large-scale cell shape transformations that is complementary to blebbing, where the plasma membrane detaches from the actin cortex and is initially unsupported when the bleb extends as a result of cytosolic pressure. Our findings provide insight into the mechanisms that drive the rapid morphological changes that occur in many physiological contexts, such as amoeboid migration and cytokinesis

Simulating Biochemical Signaling Networks in Complex Moving Geometries

Signaling networks regulate cellular responses to environmental stimuli through cascades of protein interactions. External signals can trigger cells to polarize and move in a specific direction. During migration, spatially localized activity of proteins is maintained. To investigate the effects of morphological changes on intracellular signaling, we developed a numerical scheme consisting of a cut cell finite volume spatial discretization coupled with level set methods to simulate the resulting advection-reaction-diffusion system. We then apply the method to several biochemical reaction networks in changing geometries. We found that a Turing instability can develop exclusively by cell deformations that maintain constant area. For a Turing system with a geometry-dependent single or double peak solution, simulations in a dynamically changing geometry suggest that a single peak solution is the only stable one, independent of the oscillation frequency. The method is also applied to a model of a signaling network in a migrating fibroblast

Biophysical Model of Ion Transport across Human Respiratory Epithelia Allows Quantification of Ion Permeabilities

Lung health and normal mucus clearance depend on adequate hydration of airway surfaces. Because transepithelial osmotic gradients drive water flows, sufficient hydration of the airway surface liquid depends on a balance between ion secretion and absorption by respiratory epithelia. In vitro experiments using cultures of primary human nasal epithelia and human bronchial epithelia have established many of the biophysical processes involved in airway surface liquid homeostasis. Most experimental studies, however, have focused on the apical membrane, despite the fact that ion transport across respiratory epithelia involves both cellular and paracellular pathways. In fact, the ion permeabilities of the basolateral membrane and paracellular pathway remain largely unknown. Here we use a biophysical model for water and ion transport to quantify ion permeabilities of all pathways (apical, basolateral, paracellular) in human nasal epithelia cultures using experimental (Ussing Chamber and microelectrode) data reported in the literature. We derive analytical formulas for the steady-state short-circuit current and membrane potential, which are for polarized epithelia the equivalent of the Goldman-Hodgkin-Katz equation for single isolated cells. These relations allow parameter estimation to be performed efficiently. By providing a method to quantify all the ion permeabilities of respiratory epithelia, the model may aid us in understanding the physiology that regulates normal airway surface hydration

A Whole-Body Model for Glycogen Regulation Reveals a Critical Role for Substrate Cycling in Maintaining Blood Glucose Homeostasis

Timely, and sometimes rapid, metabolic adaptation to changes in food supply is critical for survival as an organism moves from the fasted to the fed state, and vice versa. These transitions necessitate major metabolic changes to maintain energy homeostasis as the source of blood glucose moves away from ingested carbohydrates, through hepatic glycogen stores, towards gluconeogenesis. The integration of hepatic glycogen regulation with extra-hepatic energetics is a key aspect of these adaptive mechanisms. Here we use computational modeling to explore hepatic glycogen regulation under fed and fasting conditions in the context of a whole-body model. The model was validated against previous experimental results concerning glycogen phosphorylase a (active) and glycogen synthase a dynamics. The model qualitatively reproduced physiological changes that occur during transition from the fed to the fasted state. Analysis of the model reveals a critical role for the inhibition of glycogen synthase phosphatase by glycogen phosphorylase a. This negative regulation leads to high levels of glycogen synthase activity during fasting conditions, which in turn increases substrate (futile) cycling, priming the system for a rapid response once an external source of glucose is restored. This work demonstrates that a mechanistic understanding of the design principles used by metabolic control circuits to maintain homeostasis can benefit from the incorporation of mathematical descriptions of these networks into “whole-body” contextual models that mimic in vivo conditions